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3.  Statement  of  the  Problem  Studied 

3.1.  Overall  goal. 

The  overall  goal  of  the  research  project  was  to  use  kinetics,  mass  spectrometry,  and  computational 
modeling  to  investigate  mechanisms  of  aging  of  serine  hydrolases  inhibited  by  organophosphorus  (OP) 
compounds. 

3.2.  Objectives. 

The  objectives  of  the  research  were  to  test  the  following  main  hypotheses  that  were  based  on  considerable 
literature  background  and  preliminary  data: 

3.2.1.  Acetylcholinesterase  (AChE),  butyrylcholinesterase  (BChE),  and  neuropathy  target  esterase 
(NTE)  catalytic  domain  (NEST)  inhibited  by  diisopropylphosphorofluoridate  (DFP)  undergo 
postinhibitory  reactions  (e.g.,  aging)  to  yield  monoisopropylphosphyl  adducts  on  their  active  site  serines. 

3.2.2.  AChE  and  BChE  inhibited  by  A,A-diisopropylphosphorodiamidofluoridate  (mipafox,  MIP),  a 

water-soluble  analogue  of  DFP,  likewise  undergo  aging  to  yield  monoisopropylphosphyl  adducts  on  their 
active  site  serines. 

3.2.3.  In  contrast,  NEST  inhibited  by  MIP  does  not  undergo  a  similar  postinhibitory  reaction  and  retains 
a  diisopropylphosphyl  adduct  on  its  active  site  serine,  in  which  case  aging  proceeds  by  reversible 
deprotonation. 

3.2.4.  Computational  modeling  may  be  used  to  gain  insight  into  known  interactions  of  OP  compounds 
with  serine  hydrolases  and  to  make  predictions  about  currently  unknown  interactions  of  existing  or  new 
compounds  with  existing  or  new  macromolecular  targets. 

4.3.  Specific  aims. 

The  hypotheses  were  tested  and  additional  relevant  data  were  generated  by  carrying  out  the  following 
specific  aims: 

4.3.1.  Kinetics.  Bimolecular  rate  constants  of  inhibition  {k^  and  the  first-order  rate  constants  of  aging 
(^4)  were  determined  for  human  recombinant  AChE,  horse  serum  BChE,  and  human  recombinant  NEST 
inhibited  by  DFP  and  MIP.  Inhibitory  and  postinhibitory  kinetics  studies  were  needed  to  fill  gaps  in  the 
existing  database  and  especially  to  confirm  preliminary  observations  that  the  results  of  postinhibitory 
kinetics  of  MIP-inhibited  NEST  were  not  consistent  with  a  conventional  aging  mechanism. 

4.3.2.  Mass  spectrometry.  Adducts  on  the  active  site  serines  of  AChE,  BChE,  and  NEST  inactivated 
with  DFP  or  MIP  were  determined,  allowing  sufficient  time  for  aging  or  other  postinhibitory  reactions  to 
occur  as  determined  from  the  kinetics  studies  in  Aim  (1)  above.  Peptide  mass  mapping  with  surface- 
enhanced  laser  desorption/ionization  time-of-flight  mass  spectrometry  (SELDI-TOF-MS)  was  used  to 
identify  adducts. 

4.3.3.  Computational  molecular  modeling.  Computational  investigations  were  carried  out,  focusing  on 
developing  a  homology  model  of  NEST,  given  that  tertiary  structures  of  AChE  and  BChE  were  already 
known,  but  that  the  tertiary  structure  of  NTE  and  NEST  were  unknown. 
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5.  Summary  of  the  Most  Important  Results 

5.1.  Kinetics.  Rate  constants  for  inhibitory  and  postinhibitory  reactions  of  DFP  and  MIP  with  serine 
esterases  are  summarized  below  for  AChE  (Table  1),  BChE  (Table  2),  and  NEST  (Table  3): 


Table  1.  Kinetic  constants  for  inhibition  and  aging  of  AChE  by  DFP  or  MIP 

Inhibitor 

ki  (M'^min'^) 

(min'^) 

DFP 

(1.72  ±  0.07)  X  10*^ 

(2.13  ±  0.16)  X  10”^ 

MIP 

(1.47  + 0.07)  X  10^ 

(1.92  +  0.13)  X  10"^ 

Data  from  Kropp  and  Richardson  (2006). 

^  Values  are  means  ±  SEM  {n  =  A  experiments). 

Table  2.  Kinetic  constants  for  inhibition  and  aging  of  BChE  by  DFP  or  MIP 

Inhibitor 

ki  (M"'min"') 

^4  (min'^) 

DFP 

(1.83  ±  0.18)  X  10® 

(8.49  ±  0.99)  X  10"^ 

MIP 

(1.28  + 0.05)  X  10® 

(1.21+0.28)  X  10"^ 

Data  from  Kropp  and  Richardson  (2007). 

^  Values  are  means  ±  SEM  {n  =  A  experiments) 

Table  3.  Kinetic  constants  for  inhibition  and  aging  of  NEST  by  DFP  or  MIP 

Inhibitor 

ki  (M'^min'^) 

^4  (min'^) 

pH  8.0 

pH  8.0 

pH  5.2 

DFP 

(1.72  ±  0.0180)  X  lO"* 

(1.08  + 0.41)  X  10-‘ 

(1.81  +0.34)  X  lO'* 

MIP 

(1.88  ±0.061)  X  10^ 

c 

d 

""  Data  from  Kropp  et  al.  (2004). 

^  Values  are  means  ±  SEM  {n  =  A  experiments). 

^  No  activity  was  restored  for  all  time  points  (0-18  h). 

^  Activity  was  completely  restored  for  all  time  points. 

DFP  was  over  100-fold  more  potent  than  MIP  as  an  inhibitor  of  AChE  and  about  9-fold  more  potent  than 
MIP  against  NEST,  whereas  these  inhibitors  were  approximately  equipotent  against  BChE. 

Rates  of  aging  of  AChE  inhibited  by  DFP  or  MIP  were  approximately  equal,  with  half-times  of  about 
32.5  min  and  36.1  min,  respectively.  Similarly,  rates  of  aging  of  BChE  inhibited  by  DFP  or  MIP  were 
comparable,  with  half-times  of  81.3  min  and  57.3  min,  respectively.  NEST  inhibited  with  DFP  aged 
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relatively  rapidly  at  pH  8.0  or  pH  5.2  with  half-times  of  6.41  and  3.83  min,  respectively.  NEST  inhibited 
with  MIP  appeared  to  age  instantaneously  at  pH  8.0  and  no  activity  could  be  restored  during  0  to  18  h, 
whereas  at  pH  5.2  no  aging  took  place  and  all  activity  could  be  restored  for  all  time  points.  Thus,  aging  of 
MIP-inhibited  NEST  was  pH-dependent  and  the  results  were  consistent  with  aging  involving 
deprotonation  of  a  phosphoroamido  nitrogen. 

5.2.  Mass  spectrometry.  Average  masses  of  intact  (nonaged)  and  kinetically  aged  adducts  on  the 
active  site  peptides  of  serine  hydrolases  reacted  with  DFP  or  MIP  are  shown  for  AChE  (Table  4)  BChE 
(Table  5)  and  NEST  (Table  6). 


Table  4.  Differences  in  average  mass  of  MH^  ions  of  active  site  peptides  from  intact  or  aged  adducts  of 

AChE  with  DFP  or  MIP^^ 


Data  type 

Am  DFP  intact 

Am  DFP  aged 

Am  MIP  intact 

Am  MIP 
singly  aged 

Am  MIP 
doubly  aged 

Observed 

165.0  ±0.9 

123.810.7 

163.41 1.0 

ND" 

80.7  1 0.¥ 

Theoretical 

164.1 

122.1 

162.2 

121.1 

80.0 

"^Data  from  Kropp  and  Richardson  (2006). 

^Values  are  means  ±  SEM  {n>A  experiments).  The  m/z  values  (average  mass)  for  the  MH^  ion  of  the 
unreacted  tryptic  digest  peptide  containing  the  active  site  Ser  are  4270.8  ±  0.3  (observed)  and  4270.8 
(theoretical).  Theoretical  values  of  Am  (m/z)  take  into  account  the  average  mass  of  the  Ser  hydroxyl  proton 
that  is  lost  in  the  organophosphorylation  reaction.  Theoretical  nonaged  (intact)  adducts  for  DFP  and  MIP 
are  diisopropylphosphoryl  and  A,A-diisopropylphosphoroamidoyl,  respectively.  Theoretical  aged  adducts: 
monoisopropylphosphoryl  for  DFP;  singly  aged  A-monoisopropylphosphoroamidoyl  and  doubly  aged 
H2PO3  for  MIP.  Observed  m/z  and  Am  values  were  not  different  from  corresponding  theoretical  values 

(P>0.1). 

^No  peak  corresponding  to  a  singly  aged  adduct  of  MIP  was  detected. 

^  The  peak  corresponding  to  the  doubly  aged  active  site  peptide  was  confirmed  by  immunoprecipitation  by 
a  phosphopeptide  antibody  followed  by  SELDI-TOF-MS  to  obtain  a  single  peak  with  the  expected  m/z 
value  in  the  mass  spectrum. 


Table  5.  Differences  in  average  mass  of  MH^  ions  of  active  site  peptides  from  intact  or  aged  adducts  of 

BChE  with  DFP  or  MIP"’* 


Data  type 

Am  DFP  intact 

Am  DFP  aged 

Am  MIP  intact 

Am  MIP  aged 

Observed 

165.010.5 

123.410.7 

162.910.4 

122.410.7 

Theoretical 

164.1 

122.1 

162.2 

121.1 

"^Data  from  Kropp  and  Richardson  (2007). 

^Values  are  means  ±  SEM  (^  >  10  experiments).  The  m/z  values  (average  mass)  for  the  MH^  ion  of  the 
unreacted  tryptic  digest  peptide  containing  the  active  site  Ser  are  2200.5  ±0.1  (observed)  and  2200.5 
(theoretical).  Theoretical  values  of  Am  (m/z)  take  into  account  the  average  mass  of  the  Ser  hydroxyl  proton 
that  is  lost  in  the  organophosphorylation  reaction.  Theoretical  nonaged  (intact)  adducts  for  DFP  and  MIP 
are  diisopropylphosphoryl  and  A,A -diisopropylphosphoroamidoyl,  respectively.  Theoretical  aged  adducts 
for  DFP  and  MIP  are  monoisopropylphosphoryl  and  A-monoisopropylphosphoroamidoyl,  respectively. 
Observed  m/z  values  were  not  different  from  corresponding  theoretical  values  (p  >  0.1). 
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Table  6.  Differences  in  average  mass  of  MH"^  ions  of  active  site  peptides  from  intact  or  aged  adducts  of 

NEST  with  DFP  or  MIP'’  * 


Data  type 

Am  DFP  intact 

Am  DFP  aged 

Am  MIP  intact 

Am  MIP  aged’' 

Observed 

ND'^ 

123.0  ±0.6 

162.9  ±0.4 

162.9  ±0.4 

Theoretical 

164.1 

122.1 

162.2 

162.2 

‘’Data  from  Kropp  et  al.  (2004);  theoretical  values  corrected  by  subtracting  the  average  mass  for  H  (1.01 
Da). 

*  Values  are  means  ±  SEM  («  >  1 1  experiments).  The  m/z  values  (average  mass)  for  the  MIT^  ion  of  the 
unreacted  tryptic  digest  peptide  containing  the  active  site  Ser  are  2922.0  ±  0.2  (observed)  and  2922.3 
(theoretical).  Theoretical  values  of  Am  (m/z)  take  into  account  the  average  mass  of  the  Ser  hydroxyl  proton 
that  is  lost  in  the  organophosphorylation  reaction.  Theoretical  nonaged  adducts  for  DFP  and  MIP  are 
diisopropylphosphoryl  and  diisopropylphosphoroamidoyl,  respectively.  Theoretical  aged  adducts  for  DFP 
and  MIP  are  monoisopropylphosphoryl  and  A^-monoisopropylphosphoroamidoyl,  respectively.  Observed 
m/z  values  were  not  different  from  corresponding  theoretical  values  (p  >  0.1). 

^  No  peak  corresponding  to  a  classically  aged  MIP  adduct  Am  (+12 1. 1  m/v)  was  detected.  Only  a  peak 
corresponding  to  an  intact  MIP -NEST  adduct  was  detected  in  kinetically  aged  (nonreactivatable) 
MIP-NEST,  consistent  with  aged  MIP -NEST  being  the  deprotonated  species,  which  would  be  protonated 
under  the  acidic  conditions  used  for  SELDI-TOF  MS. 

^  No  peak  corresponding  to  a  nonaged  DFP  adduct  was  detected,  consistent  with  relatively  rapid  aging  of 
the  DPF-NEST  conjugate  (half-time  of  aging,  6.4  min). 


Table  4  shows  that  DFP-inhibited  AChE  aged  in  the  expected  manner  to  give  a  peak  in  the  SELDI-TOF 
MS  spectrum  corresponding  to  a  monoisopropylphosphoryl  adduct  on  the  active  site  Ser  peptide. 
However,  no  peak  corresponding  to  a  A-monoisopropylphosphoroamidoyl  adduct  on  the  active  site 
peptide  was  detected  in  samples  from  kinetically  aged  MIP-inhibited  AChE.  Instead,  a  peak 
corresponding  to  a  simple  phosphate  adduct  (H2PO3)  was  detected.  This  surprising  result  was  confirmed 
by  immunoprecipitation  using  a  phosphopeptide  antibody  followed  by  SELDI-TOF  MS.  The  resulting 
spectrum  consisted  of  a  single  peak  corresponding  to  the  active  site  peptide  with  a  phosphate  (H2PO3) 
adduct  (bracketed  by  two  peaks  from  internal  molecular  weight  standards  that  were  added  to  the  MS 
chip).  Thus,  MIP-inhibited  AChE  ages  by  loss  of  both  isopropylamine  groups,  as  depicted  in  Fig.  1. 

BChE  inhibited  by  DFP  or  MIP  aged  in  the  conventional  manner  as  shown  by  the  MS  data  in  Table  5.  In 
each  case,  the  mass  spectra  of  kinetically  aged  samples  showed  disappearance  of  the  peak  corresponding 
to  the  intact  OP  adduct  and  appearance  of  a  peak  corresponding  to  the  monoisopropylphosphyl  adduct. 

Aging  of  DFP-inhibited  NEST  resulted  in  a  peak  corresponding  to  the  expected 
monoisopropylphosphoryl  adduct,  whereas  the  mass  spectrum  of  kinetically  aged  MIP-inhibited  NEST 
showed  only  the  intact  adduct.  This  result,  coupled  with  the  kinetics  results  showing  complete  aging  at 
pH  8.0  and  complete  reactivation  at  pH  5.2  indicate  that  aging  of  MIP-inhibited  NEST  occurs  by 
deprotonation  of  the  acidic  phosphoroamido  hydrogen. 

The  different  aging  pathways  for  MIP-inhibited  esterases  are  shown  below  in  Fig.  1.  It  is  noteworthy  that 
the  three  MIP-inhibited  esterases  underwent  aging  by  three  distinct  mechanisms  (the  details  of  which 
require  further  elucidation)  to  yield  three  different  adducts.  These  results  indicate  that  further  research  is 
needed  to  understand  more  fully  the  mechanisms  of  aging  of  phosphoroamidated  serine  esterases.  In 
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addition,  the  fact  that  phosphoroamidates  can  inhibit  esterases  to  yield  different  adducts  has  important 
implications  for  the  delineation  of  biomarkers  of  exposure  to  OP  compounds. 
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Figure  1.  Pathways  of  inhibition  and  aging  of  MIP  for  three  serine  esterases.  (A)  NEST,  (B)  AChE,  and 
(C)  BChE.  Inhibition  of  each  serine  hydrolase  occurs  by  nucleophilic  attack  of  the  Ser  hydroxyl  on  the 
phosphorus  of  MIP,  displacing  the  primary  leaving  group,  fluoride,  and  yielding  an  A,A^-diisopropyl- 
phosphorodiamido  adduct,  which  then  undergoes  aging  by  a  different  mechanism  for  each  enzyme. 
Although  the  exact  mechanisms  (SnI  or  Sn2)  of  aging  for  MIP-inhibited  AChE  and  BChE  require 
definitive  elucidation,  results  to  date  are  consistent  with  dual  (AChE)  or  single  (BChE)  Sn2  displacement 
of  isopropylamine  by  water  as  shown  here.  The  charge  of  each  aged  moiety  is  shown  as  expected  for  pH 
8.0,  although  these  species  would  be  expected  to  be  protonated  under  the  acidic  conditions  used  for 
SELDI-TOF  MS.  (Reproduced  from  Kropp  and  Richardson,  2007). 

5.3.  Computational  molecular  modeling. 

5.3.1.  Docking  studies  of  isomalathion  stereoisomers  and  Torpedo  AChE.  Early  computational 
molecular  modeling  work  was  carried  out  to  investigate  stereoselective  inhibition  of  Torpedo  AChE  and 
other  serine  esterases  by  isomalathion  stereoisomers  as  illustrated  in  Fig.  2.  The  results  of  computational 
docking  studies  were  consistent  with  the  experimental  results  obtained  from  kinetics  and  SELDI-TOF 
MS.  Namely,  that  diethyl  thiosuccinyl  was  the  leaving  group  for  (li?)-isomers  and  thiomethyl  was  the 
leaving  group  for  (15)-isomers  of  isomalathion  (Doom  et  al.,  2003) 

5.3.2.  Homology  model  of  the  patatin  domain  ofNTE.  The  three-dimensional  stmcture  of  NTE  is 
unknown.  However,  bioinformatics  studies  revealed  that  the  catalytic  (esterase)  domain  had  sequence 
homology  to  patatinl?  (patl7),  a  lipase  found  in  potatoes  and  other  plants  whose  crystal  stmcture  had 
been  determined  (Rydel  et  al.,  2003).  Using  the  stmcture  of  patl7  as  a  template,  we  constmcted  a 
homology  model  of  the  esterase  domain  ofNTE  (PNTE)  as  shown  in  Fig.  3  (Wijeyesakere  et  al.,  2007; 
Wijeyesakere,  2009). 
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Figure  2.  Results  of  computer-based  docking  experiments  for  isomalathion  enantiomers  and  Torpedo 
AChE.  (A)  (IR,3R)-  and  (B)  (15',3A)-isomalathion  in  the  active  site  of  Torpedo  AChE.  The  isomalathion 
enantiomers  are  in  ball-and-stick  format,  whereas  Torpedo  AChE  residues  are  displayed  as  sticks.  Colors 
were  assigned  to  atoms  as  follows:  C,  gray;  H,  green;  N,  blue;  O,  red;  P,  purple;  and  S,  yellow.  Only  the 
catalytic  triad  and  active  site  residues  in  close  proximity  to  isomalathion  are  shown.  The  catalytic  Ser  200 
is  located  directly  beneath  the  isomalathion  phosphorus  atom.  Diethyl  thiosuccinyl  is  predicted  to  be 
positioned  apically  to  the  active  site  Ser  200  for  (li?,3i?)-isomalathion.  In  contrast,  thiomethyl  is  predicted 
to  be  positioned  apically  to  Ser  200  for  the  (liS,3A)-enantiomer.  The  active  site  entrance  is  located  directly 
above  the  isomalathion  molecule  in  each  panel  as  indicated  by  the  arrow  (Reproduced  from  Doom  et  al., 
2003). 
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Figure  3.  (A)  Sequence  alignment  of  PNTE  and  pat  17.  Identical  residues  are  boxed  in  black;  homologous 
residues  in  grey.  Secondary  structures  of  PNTE  and  pat  17  are  shown  along  with  their  respective  primary 
structures  as  arrows  (P-strands)  and  rectangles  (a-helices).  a-helices  and  P-strands  in  the  PNTE  sequence 
are  numbered  al-a7  and  B1-B6,  respectively.  (B)  Ribbon  diagram  of  the  putative  tertiary  structure  for 
PNTE.  Blue  arrow  denotes  the  location  of  the  nucleophilic  elbow.  (C)  Organization  of  the  predicted 
active  site  of  PNTE,  showing  the  proposed  Ser966-Aspl086  catalytic  dyad.  Ser966,  Asp960,  and 
Asp  1086  are  depicted  as  capped  sticks.  Ser966  OY-Asp960  and  052  Ser966  Oy-Asp  1086  052  distances 
are  shown  as  dotted  lines.  (Reproduced  from  Wijeyesakere  et  aL,  2007  with  corrections  in  the  sequence 
numbering  from  Wijeyesakere,  2009). 
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A  major  finding  from  the  homology  model  of  PNTE  was  that  the  active  site  was  predicted  to  be  a  Ser966- 
Aspl086  catalytic  dyad  characteristic  of  certain  lipases  (Wijeyesakere  et  al.,  2007)  rather  than  the 
previously  proposed  Ser966-Asp960-Aspl086  catalytic  triad  (Atkins  and  Glynn,  2000). 

5.3.3  Docking,  simulated  inhibition,  and  simulated  aging  using  the  PNTE  model. 

Docking  studies  were  undertaken  with  DFP  and  MIP,  using  the  PNTE  model  as  the  target.  The  structures 
of  DFP  and  MIP  were  generated  in  Chem3D  (CambridgeSoft,  Cambridge,  MA)  and  minimized  using  the 
MM2  protocol  to  a  root  mean  square  deviation  (RMSD)  gradient  of  0.01  prior  to  their  use  in  docking. 
Docking  of  DFP  into  the  active  site  of  PNTE  was  undertaken  using  the  automated  docking  algorithm  in 
Fujitsu  CAChe  (Fujitsu  BioSciences,  Sunnyvale,  CA)  employing  a  potential-mean-force  technique  with  a 
Lamarckian  genetic  algorithm  (run  over  3,000  generations  with  a  population  size  =  50,  crossover  rate  = 

0.8  and  mutation  rate  =  0.3)  (Morris  et  al.,  1998).  Flexibility  was  permitted  in  the  ligands.  The  resulting 
models  were  analyzed  in  CAChe  and  PyMol  (DeLano,  2002).  Similar  techniques  were  undertaken  to  dock 
MIP  into  the  active  site  of  PNTE  in  order  to  compare  its  binding  affinity  to  that  of  DFP.  This  technique 
was  similar  to  that  carried  out  in  our  laboratory  on  other  serine  hydrolases  and  stereoisomers  of 
isomalathion  (although  in  those  cases  flexibility  was  allowed  in  both  the  ligand  and  a  region  of  the 
macromolecule)  (Doom  et  al.,  2001,  2003). 

As  shown  in  Fig.  4,  Docking  of  DFP  and  MIP  into  the  active  site  of  PNTE  revealed  that  the  phosphoms 
atom  of  DFP  was  located  4.7  A  from  Ser966  OG  and  3.7  A  from  Asp  1086  OD2,  while  the  phosphoms 
atom  of  MIP  was  located  5.2  A  from  Ser966  OG  and  4.2  A  from  Aspl086  OD2.  These  results  indicated 
that  these  OP  compounds  could  be  accommodated  in  the  active  site  of  the  PNTE  model  in  an  orientation 
suitable  for  organophosphorylation.  Moreover,  the  binding  energies  of  DFP  and  MIP  to  the  active  site  of 
PNTE  were  calculated  to  be  -36.1  kcal/mol  and  -30.77  kcal/mol,  respectively,  thereby  indicating 
favorable  energetics  for  interactions  between  these  ligands  and  the  active  site  of  PNTE  (Wijeyesakere, 
2009). 


Figure  4.  Docking  DFP  and  MIP  into  the  PNTE  model.  Docking  of  (A)  DFP  and  (B)  MIP  into  the  active 
site  of  PNTE.  The  phosphoryl  oxygens  of  the  OP  compounds  make  polar  contact  with  Arg941  in  the 
oxyanion  hole  with  the  Arg941  NHl  being  3.0  A  and  3.2  A  from  the  phosphoryl  oxygens  of  DFP  and 
MIP,  respectively.  (Reproduced  from  Wijeyesakere,  2009). 

Analysis  of  the  primary  sequence  of  PNTE  had  previously  indicated  that  Arg941  would  be  a  potential 
candidate  residue  in  the  oxyanion  hole  formed  by  Gly938,  Gly939,  Ala940,  and  Arg941  (Wijeyesakere  et 
al.,  2007).  Indeed,  the  results  of  the  docking  studies  revealed  that  the  NHl  of  Arg941was  located  3.0  A 
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and  3.2  A  from  the  phosphoryl  oxygen  of  DFP  and  MIP,  respectively  (Fig.  4).  These  results  lend  support 
to  a  role  of  Arg941  in  the  PNTE  oxyanion  hole  and  suggest  that  this  residue  stabilizes  the  OP  ligand  as  it 
enters  the  PNTE  active  site  in  preparation  for  nucleophilic  attack  on  the  OP  phosphorus  by  the  Ser966 
OG  (Wijeyesakere,  2009). 

Simulated  inhibition  and  simulated  aging  studies  were  undertaken  with  aged  and  non-aged  DFP  (MIP  was 
excluded  in  this  portion  of  the  study  due  to  its  unusual  mechanism  of  aging  via  deprotonation  as  opposed 
to  the  loss  an  alkyl  side-chain  (Kropp  et  al.,  2004).  This  allowed  for  testing  in  silico  the  hypothesis  that 
aging  results  in  a  conformational  change  in  PNTE  that  is  distinct  from  the  conformational  change  induced 
solely  by  inhibition.  A  similar  computational  procedure  has  been  carried  out  recently  through  virtual 
mutation  of  a  serine  to  phosphoserine  (Groban  et  al.,  2006). 

Simulated  inhibition  was  undertaken  using  the  structure  of  diisopropylphosphoryl-serine  (DPSer)  that  was 
built  from  the  structure  of  monoisopropylphosphoryl-serine  (MPSer)  downloaded  from  the  Hetero¬ 
compound  Information  Centre  (http://alpha2.bmc.uu.se/hicup/;  Kleywegt  and  Jones,  1998)  using 
ChemSD.  Ser966  in  PNTE  was  replaced  with  DPSer  using  O  10.1  (Jones  et  aL,  1991).  The  lego-loop 
procedure  in  O  was  employed  to  fix  chain  breaks  arising  from  this  substitution.  Hydrogens  were  added  to 
the  DPSer-PNTE  model  at  pH  7.0  followed  by  300  rounds  of  minimization  in  Fujitsu  CAChe  using  the 
MM3  protocol.  Simulated  aging  of  DPSer-PNTE  was  undertaken  by  deleting  an  isopropyl  group  from 
DPSer,  thereby  adding  a  negative  charge  to  the  resulting  MPSer  adduct.  The  MPSer-PNTE  structure  was 
minimized  in  Fujitsu  CAChe  using  the  MM3  protocol  as  described  for  simulated  inhibition.  The  resulting 
structures  were  analyzed  in  CAChe  and  PyMol.  Packing  of  residues  around  the  PNTE  active  site  in  the 
minimized  DPSer-PNTE  and  MPSer-PNTE  models  was  undertaken  via  visual  assessment. 

The  crystallography  and  NMR  system  (CNS;  Brunger  et  al,  1998)  ver  1.2  was  utilized  to  calculate  the 
root-mean  square  deviation  (RMSD)  between  the  Ca  backbones  of  native  PNTE  and  DPSer-PNTE  as 
well  as  the  RMSD  between  the  backbones  of  native  PNTE  and  MPSer-PNTE.  The  calculated  RMSD 
values  were  employed  as  a  measure  of  the  degree  of  overall  conformational  change  arising  from  SI  and 
SA  of  NTE  by  a  neuropathic  OP. 

Fig.  5  A  shows  that  simulated  inhibition  resulted  in  a  predicted  expansion  around  the  active  site  of  PNTE, 
with  the  distance  between  the  Ser966  OG  and  Aspl086  OD2  increasing  to  5.2  A,  compared  to  2.9  A  in 
the  native  PNTE  structure.  Simulated  aging  of  the  covalently  bound  DFP  adduct  resulted  in  a  slight 
additional  expansion  around  the  active  site  relative  to  the  inhibited  state  with  the  Ser966  OG  and 
Aspl086  OD2  located  5.5  A  apart  (compared  to  the  2.9  A  distance  seen  in  the  native  structure)  (Fig.  5B). 
This  rotamer  change  is  predumably  due  to  charge  repulsion  betweent  the  negatively  charged  MPSer 
moiety  and  the  negatively  charged  aspartate. 

The  negatively  charged  oxygen  in  the  aged  DFP  adduct  makes  polar  contact  with  Arg941  NHl  (2.7  A), 
Gly938  N  (3.4  A)  and  Gly939  N  (2.8  A),  which  form  the  oxyanion  hole  of  PNTE  (Fig.  6).  Thus,  it  can  be 
seen  that  the  negative  charge  on  the  aged  OP  adduct  is  stabilized  by  the  oxyanion  hole  in  a  manner  similar 
to  that  seen  in  normal  catalysis  with  a  carboxylic  acid  ester  tetrahedral  intermediate.  This  is  also  similar  to 
the  manner  in  which  the  oxyanion  hole  of  AChE  aids  in  the  stabilization  of  the  aged  DFP  adduct  in  the 
crystal  structure  of  AChE  (PDB  ID  2dfp). 

It  was  evident  that  covalent  binding  of  DFP  to  Ser966  residue  did  not  result  in  any  significant 
conformational  change  in  the  overall  tertiary  structure  of  PNTE  (Fig.  7A).  The  calculated  Ca  backbone 
deviation  of  native  PNTE  versus  its  DFP-inhibited  state  was  1 .4  A.  In  addition,  simulated  aging  of  the 
inhibited  PNTE  structure  did  not  show  any  significant  conformational  change  relative  to  its  native  and 
inhibited  states  (Fig.  7B).  The  calculated  backbone  RMSD  between  native  PNTE  versus  its  aged  state  was 
1.6  A.  The  backbone  Ca  deviation  between  the  inhibited  and  aged  states  of  PNTE  was  0.2  A. 
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Figure  5.  PNTE  active  site  after  simulated  inhibition  and  simulated  aging  with  DFP.  (A).  The  model 
predicts  that  inhibition  results  in  the  active  site  Ser966  OG  moving  5.1  A  away  from  Asp  1086  OD2 
(compared  to  2.9  A  between  these  residues  in  the  native  PNTE  structure).  (B).  Aging  is  predicted  to  result 
in  a  slight  additional  change,  with  the  active  site  residues  moving  5.5  A  apart  (presumably  due  to  charge 
repulsion  by  the  negatively  charged  monoisopropylphosphoryl-Ser966  and  Asp  1086).  However,  these 
changes  represent  different  rotamer  conformations  rather  than  perturbations  in  the  amino  acid  backbones 
of  these  residues.  (Reproduced  from  Wijeyesakere,  2009). 


Arg940 


Ser966 


Figure  6.  Simulated  stabilization  of  the  aged  DFP-PNTE  conjugate.  The  negative  charge  on  the  DFP- 
Ser966  adduct  is  stabilized  by  residues  of  the  oxyanion  hole  (Gly938,  Gly939  and  Arg940),  in  a  manner 
similar  to  what  is  expected  for  stabilization  of  the  tetrahedral  intermediate  of  an  ester  during  normal 
catalysis  and  similar  to  that  seen  the  crystal  structures  of  other  aged  OP-serine  esterase  conjugates.  The 
negatively  charged  oxygen  is  indicated  along  with  contact  distances  between  it  and  Gly938  N  (3.4A), 
Gly939  N  (2.8  A)  and  Arg940  NHl  (2.7  A).  (Reproduced  from  Wijeyesakere,  2009). 
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Figure  7.  Simulated  structural  changes  of  PNTE  Ca  backbone  after  inhibition  and  aging.  The  Ca 
backbone  structures  of  native  PNTE  (green)  is  compared  to  the  backbone  structures  of  PNTE  following 
(A)  inhibition  by  DFP  (purple),  or  (B)  aging  of  the  DFP  adduct  (orange).  Inhibition  of  PNTE  does  not 
result  in  a  significant  overall  conformational  change  relative  to  its  native  state  (Ca  backbone  RMSD  = 
1.4  A).  Similarly,  aging  is  not  associated  with  a  significant  conformational  change  relative  to  native 
PNTE  (Ca  backbone  RMSD  =1.6  A).  (Reproduced  from  Wijeyesakere  et  aL,  2009). 

These  findings  can  be  compared  to  previous  molecular  dynamics  research  on  AChE  undertaken  by  Hurley 
et  al.  (2005)  that  reported  a  net-RMSD  of  1.2  A  associated  with  in  silico  inhibition  (but  not  aging)  of 
AChE  by  the  nerve  agent  VX.  This  RMSD  is  comparable  to  our  result  of  1.4  A  for  simulated  inhibition  of 
PNTE  by  DFP.  With  respect  to  aging  of  an  OP  adduct,  the  X-ray  crystal  structures  of  Torpedo  californica 
AChE  inhibited  and  aged  by  DFP,  sarin,  and  soman  (PDB  IDs  2dfp,  Icfj,  and  Isom,  respectively) 
reported  no  significant  conformational  changes  in  the  overall  structure  of  the  enzyme,  although  some 
localized  conformational  changes  were  reported  in  the  AChE  acyl  pocket  (Millard  et  al.,  1999).  Similarly, 
analysis  of  the  crystal  structure  of  native  human  BChE  (PDB  ID  IpOi)  with  its  DFP-  and  echothiophate- 
aged  forms  (PDB  ID  Ixlu  and  Ixlv,  respectively)  shows  no  evidence  for  a  global  conformational  change 
associated  with  aging  (Nachon  et  al.,  2005).  The  Ca  RMSDs  for  the  DFP-  and  echothiophate-aged  forms 
of  BChE  relative  to  the  native  form  of  this  enzyme  are  0.51  A  (DFP-aged  vs.  native)  and  0.55  A 
(echothiophate-aged  vs.  native).  This  lack  of  a  global  conformational  change  following  aging  of  an  OP- 
esterase  conjugate  is  consistent  with  the  findings  of  our  molecular  modeling  effort  of  the  catalytic  domain 
of  NTE  and  serves  to  suggest  that  aging  of  NTE  is  not  expected  to  result  in  a  significant  change  in  the 
structure  of  the  enzyme  relative  to  its  native  state. 

Thus,  if  aging  of  NTE  results  in  a  toxic  gain  of  function  of  the  protein  that  is  responsible  for  initiation  of 
neurodegeneration  (Glynn,  2000),  then  our  modeling  results  on  PNTE  indicate  that  such  a  change  does 
not  arise  from  a  global  conformational  change.  On  the  other  hand,  initiation  of  neurodegeneration  could 
arise  from  a  loss  of  function  of  NTE,  such  as  its  proposed  role  in  maintaining  the  homeostasis  of 
membrane  phospholipids  (van  Tienhoven  et  al.,  2002;  Vose  et  al.,  2008;  Zaccheo  et  al.,  2004).  Further 
research  involving  modeling  and  experimental  studies  on  full-length  NTE  and  related  serine  hydrolases 
will  be  required  to  address  these  important  questions. 
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